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^So neutron pairing vs observations of cooling neutron stars 
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Abstract. As shown recently by Kaminker et al. (2001), current observations of thermal emission of isolated 



middle-aged neutron stars (NSs) can be explained by cooling of NSs of different masses with the cores composed 
of neutrons, protons and electrons, assuming rather strong superfluidity (SF) of protons, weak triplet-state SF 
of neutrons and neglecting singlet-state SF of neutrons. We show that this explanation remains correct in the 
presence of singlet-state SF of neutrons in the NS crust and outermost core but under stringent constraints on 
the density profile of the SF critical temperature rcns5(p). In order to explain observations of (young and hot) RX 
J0822-43 and (old and warm) PSR 1055-52 and RX J1856-3754 as cooling not too massive NSs, the maximum 
Tens'' should be rather high (^ 5 x 10^ K) and the decrease of Tcus{p) outside the maximum should be sharp. 
These results place important constraints on the models of nucleon SF in dense matter. 
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1. Introduction 

Cooling of NSs depends on the poorly known properties of 
matter in NS interiors. Combined with the best observa- 
tional data, the cooling theory can be used to explore the 
properties of this matter, in particular, critical tempera- 
tures Tc of SF of neutrons (n) and protons (p) (supple- 
menting very model-dependent microscopic calculations, 
e.g., Lombardo & Schulze 2001 ). SF of nucleons reduces 
the emissivity of neutrino reactions and affects the nu- 
cleon heat capacity (e.g., Yakovlev et al 
SF initiates a specific neutrino emission associated with 
Cooper pairing of nucleons (Flowers et al. 1976). In this 



plained with the models of cooling NSs of different masses 
assuming rather strong proton SF and weak neutron 
SF. The models require neither exotic composition of NS 
core nor additional NS reheating. 

We extend the results of Paper I by including ^Sq neu- 
tron superfluidity and show that, contrary to our expec- 
tation, this superfluidity is important. 



1999| ). Moreover, 2. Cooling models with ^Sq neutron pairing 



way nucleon SF becomes a strong regulator of NS cooling. 

Recently Kaminker et al. ( 2001 , hereafter Paper I) pro- 
posed the interpretation of observations of thermal emis- 
sion from eight isolated middle-aged NSs (Fig. 1) using 
simple models of cooling NSs with superfluid cores (com- 
posed of neutrons, protons and electrons). Accordingly, 
one can consider three types of nucleon SFs with density 
dependent Tc{p): singlet-state (^Sq) neutron SF in the NS 
inner crust and the outermost core (Tc — Tens); triplet- 
state (^P2) neutron SF in the core (Tent); and ^Sq proton 
SF in the core (Tcp). For simplicity, in Paper I ^Sq neutron 
SF was neglected. The main reason was that, being mainly 
located in a thin crust, this SF seemed unimportant for 
middle-aged NSs. In Paper I the observations were ex- 
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We simulate NS cooling with our fully relativistic non- 
isothermal cooling code (e.g., Gnedin et al. 2001). The 
physics input is mostly the same as in Paper I. In par- 
ticular, we adopt a moderately stiff equation of state 
(EOS) of matter in the NS core proposed by Prakash et 
al. ( |1988 ) (their model I with the compression modulus 
of saturated nuclear matter K = 240 MeV). The outer 
heat blanketing NS envelope is assumed to be composed 
of iron, but the atmosphere may be made of hydrogen. 
The maximum mass of our NS models is 1.977 Mq (cen- 
tral density /Oc,max = 2.575 x 10^^ g cm~'^). For the given 
EOS, a powerf ul dire ct Urea process of neutrino emission 
(Lattimer et al. |l99l| ) is forbidden if M < Mu = 1.358 M© 
{pc < 7.851 X 10^^ g cm''^). In this case, a NS undergoes a 
slow cooling mainly via modified Urea process. A NS with 
M > Md possesses a central kernel, where direct Urea 
process is open, initiating fast cooling. 
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Fig. 1. Observational limits on surface temperatures of 
eight NSs compared with cooling curves for NS models 
with masses from 1.35 to 1.5 Mq. Dot-and-dash curves 
are obtained including proton SF in the NS core shown in 
Fig. 2. Solid curves include, in addition, model 1 of Sq 
neutron SF from Fig. 2. 



We will confront theoretical cooling curves with the 
same observational data as in Paper I; the only exclu- 
sion is RX J1856-3754, for which new data have become 
available. Figure 1 shows the limits on the effective sur- 
face temperature, Tg°°, as measured by a distant observer, 
versus stellar age t for eight NSs. The youngest three are 
radio-quiet NSs in supernova remnants; the oldest is also 
a radio-quiet NS; the others are observed as radio pulsars. 
The values of are taken from t he following sources: RX 
J0822- 43 — ZavHn et al. ( |1999D , IE 1207-52 — Zavli n 
et al. ( |199^ ), RX J0002-)-62 — Zav lin fc Pavlov ( |1999D , 
PSR 0833-45 (Vela) — P avlov et al. ( ^OOlD , PSR 0656-M4 
— Possenti et al (|l996| ), PSR 0630-^178 (Gemin ga) - - 
Halpern & Wang ( |l997l ), PSR 1055-52 — O gelman ( |1995D , 
and RX J1856-3754 — Pons et al. ( ^00 ij ). As in Paper 
I the values of for the four youngest NSs are in- 
ferred from the observations using hydrogen atmosphere 
models while for other SNs (but RX J1856) they are in- 
ferred using the blackbody model. For RX J 1856, we take 
= 0.52 ± 0.07 MK, as inferred using the analytic fit 
with the Si-ash atmosphere model, and t = 9 x 10^ K 



(Pons et al. 2001). We assume that wide errorbar of 
given by this model is in line with current poor under- 
standing of t herma l emission from RX J1856 (als o wee 
Burwitz et al 



2001 



It is also possible (Pons et al. |2001 



Burwitz et al. 2001 ) that the source has a colder surface 
(0.25 MK) with a hot spot. This case would be easier for 
our interpretation (see below) than the case we adopt. 
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Fig. 2. Density dependence of the critical temperatures 
for one model of the proton SF (dots-and-dashes) and 
three models 1, 2 and 3 of ^So neutron SF (Table 1) used 
in cooling simulations (Figs. 1 and 3). Vertical dotted lines 
indicate neutron drip, core-crust interface, and the direct 
Urea threshold. 



We will not rely on any particular (very model depen- 
dent) microscopic calculation of SF critical temperature 
Tc{p). Instead, following Paper I, we parameterize Tc as 



fc2 



k^ + kl {k~k2Y + kl 



(1) 



for k < k2, and Tc ~ for k > k2- Here, k = fcpN is 
the nucleon (N=n or p) Fermi wavenumber (measured in 
fm^^). The parameters Tq, fci, fc2, ^3 of our SF models are 
given in Table 1. 

As in Paper I we adopt a typical proton SF (Fig. 2). We 
also assume that neutron SF in the NS core is weak 
(maximum T™^'' < 10* K) and does not affect cooling 
(see Paper I for details). Neglecting ^Sq neutron SF, we 
obtain the dot-and-dashed cooling curves displayed in Fig. 
1 for several values of M . The upper curve is calculated 
for M = 1.35 Mq but it is actually the same for all M 
from 1.1 M0 to Md. The dot-and-dashed curves in Fig. 

1 are slightly different from analogous curves in Paper I 
since now we use a denser grid of mass zones in the cooling 
code. In agreement with Paper I we see that tuning M we 
can explain the observations of all eight sources. 

Now let us include -^Sq neutron SF. Figure 2 shows 
the three models we consider. Models 1 (solid line) and 

2 (dashes) correspond to about the same, rather strong 
SF (with maximum Tens « 7 x 10^ K). The main differ- 
ence is that Tcns(p) in model 2 has flatter maximum and 
sharper fall in the maximum wings. Model 3 (dots) de- 
scribes much weaker SF, with maximum Tens ~ 2.4 x 10^ 
K and a narrower density proflle. 

The solid lines in Fig. 1 show cooling curves calcu- 
lated under the same assumptions as the dot-and-dashed 
lines but including ^Sq neutron SF (model 1). As expected, 
the presence of this SF shortens the initial thermal relax- 
ation stage from 100 yrs to 30-50 yrs (Gnedin et al. 2001) 
but does not affect strongly the cooling of middle-aged 
NSs with open direct Urea process (M > Mq). Thus, 
^So neutron SF does not change the results of Paper I 
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Table 1. Parameters of ^So SF models in Eq. (1) 



SF 


To/lO'-* K 


fci, fm ^ 


A;2, fm ^ 


ks, fm 


P 


20.29 


1.117 


1.241 


0.1473 


In 


10.2 


0.6 


1.45 


0.1 


2n 


7.9 


0.3 


1.45 


0.01 


3n 


1800 


21 


1.45 


0.4125 



for M > A/d and the proposed interpretation of obser- 
vations of five sources, IE 1207-52, RX J0002+62, Vela, 
PSR 0656+14, and Geminga. 

However, this SF drastically affects the cooling of low- 
mass {M < Md) NSs. Its main effects are to lower at 
t ^ 3 X 10^ yrs by an additional neutrino emission associ- 
ated with the ^So Cooper pairing of neutrons and to lower 
Tg°° further at i ^ 3 x 10^ yrs by reducing heat capacity 
of neutrons in the crust. The importance of the Cooper 
emission can be explained as follows. The neutrino emis- 
sion from the core of a middle-aged NS is strongly reduced 
by SF suppression of the majority of neutrino reactions 
(T <C Tcp). The neutrino emission due to ^Sg neutron pair- 
ing is also greatly suppressed in the bulk of the inner NS 
crust, where T <C Tcnsj but remains active in two thin lay- 
ers, where T^nsip) closer to the internal temperature T: 
near the neutron drip point and in the outermost part of 
the core. Consider, for instance, the outermost core. Since 
the emission layer is thin, the emissivity can be written 
as Q^^ = QqT"^ F{T /Tens), where Qo may be regarded as 
constant (easily obtained from Eq. 79 in Yakovlev et al. 



1999) and F{t) is a known function. If we approximate 
ens on radial coordinate r in our thin 
(r) = T + A{r - ro) 



the dependence of Tc, 

layer by a linear function Tcns('') = T + A(r — rp) and 
integrate over this layer, we obtain the Cooper-pairing 
(non-redshifted) neutrino luminosity: 



6.63, 



(2) 



where h = T/A is the characteristic width of the layer 
(variation scale of T^nsi''') at ~ ^o)- Adopting a typi- 
cal value ro = 12 km and the neutron number density 
n„ = 0.1 fm~^ we have « 8.9 x lO^^ (/i/lOO m) erg 
s"\ where Tg = T/IO^ K. For a typical scale h - 100 m, 
the luminosity V^^^ is comparable to the standard (modi- 
fied Urea) luminosity of the entire non-superfluid NS core! 
It may easily outweight the SF suppression of the neu- 
trino luminosity in other reactions. In Paper I we analyzed 
analogous effect with regard to '^P2 neutron SF in the core 
and showed that such SF, if available, strongly accelerates 
cooling (even for M > Mq) and complicates interpreta- 
tion of observations. Here we demonstrate that a similar 
but less pronounced effect is produced by ^Sq neutron SF 
at M < Md. 

According to Fig. 1, the ^Sq neutron SF compHcates 
interpretation (Paper I) of observations of the youngest 
and hot source RX J0822-43 and of two oldest and warm 
sources, PSR 1055-52 and RX J1856-3754 by the NS mod- 
els with M < Md . The less stringent complication for RX 
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Fig. 3. Cooling curves of the 1.35 Mq NS model vs ob- 
servations of RX J0822-43, PSR 1055-52 and RX J1856- 
3754. Dot-and-dashed curve: proton SF in the NS core. 
Solid, dashed and dotted curves include, in addition, mod- 
els 1, 2 and 3 of ^ So neutron SF (Fig. 2), respectively. Full 
solid line is the same as in Fig. 1. Thin solid curve acc is 
calculated assuming the presence of 1O~^°M0 of hydro- 
gen on the NS surface. Thin solid curve mag is obtained 
assuming the dipole surface magnetic field (10^^ G at the 
magnetic pole). 



J1856-3754 is evidently associated with too wide error- 
bar of Tg°°. Let us focus (Fig. 3) on cooHng of low-mass 
NS models with regard to interpretation of these three 
sources. For certainty, we again exploit the 1.35 M© model 
although nearly the same is true for any M from 1.1 M© 
to Md- To save the proposed interpretation we must rise 
the upper cooling curve in Fig. 1. 

The obvious solution is to take ""^So neutron SF with 
steeper slopes of Tcns(p) near the crust-core interface and 
the neutron drip point. This decreases the characteristic 
scale h and the "harmful" luminosity L'^^ . For example, 
taking model 2 of ^ So neutron SF instead of model 1 (Fig. 
2) wc obtain the dashed cooling curve (Fig. 3) which is 
closer to the dot-and-dashed curve than the thick solid 
curve. (Another example: shifting Tcns(p) for model 2 into 
the crust would additionally rise the cooling curves.) Note 
that the cooling curves are insensitive to the details of 
the Tcns(p) profiles near the maximum of Tcns(p) as long 
as ^ 5 X 10^ K, but they are very sensitive to 

the decreasing slopes of Tcns(p)- A similar effect was re- 
ported in Paper I with regard to the Tcp(p) profiles. Taking 
smoother and lower Tcns(p), model 3, we obtain a NS that 
is much colder than required by observations (dotted line 
in Fig. 3). Therefore, ^So neutron SF with T^f < 5 x 10^ 
K and/or with smoothly decreasing slopes of the Tcns{p) 
profile near the crust-core interface and neutron drip point 
violates the proposed interpretation. 
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Let us stress that the observations of these three 
sources can be fit even with our initial model 1 of ^Sq 
neutron SF. The high surface temperature of RX J0822- 
43 can be explained assuming the presence of low-mass 
surface envelope of light elements. This effect is modeled 
using the results of Potekhin et al. (1997). Light elements 
increase the electron thermal conductivity of NS surface 
layers and raise at a neutrino cooling stage, t ^ 3 x 10^ 
yrs (curve acc in Fig. 3). In order to explain the observa- 
tions of PSR 1055-52 and RX J1856-3754, we can as- 
sume again model 1 SF, iron surface and the dipole sur- 
face magnetic field (~ 10^^ G at the magnetic pole; line 
mag). Such a field makes the NS surface layer s overall less 
heat-transparent (Potekhin & Yakovlev 2001 ), rising 
at i ^ 3 X 10^ yrs. Note that the dipole field ^ 10" G has 
the opposite effect, similar to that produced by the surface 
envelope of light elements. Thus, we can additionally vary 
cooling curves by assuming the presence of light elements 
and/or the magnetic field on the NS surface. However, 
these variations are less pronounced than those produced 
by nucleon SF. For instance, we cannot reconcile the cool- 
ing curves with observations of PSR 1055-52 assuming 
model 3 of neutron SF with any surface magnetic field. 

3. Summary 

Extending the results of Paper I, we have shown that ob- 
servations of thermal emission from eight isolated NSs are 
consistent with a simple model of cooling superfluid NSs of 
different masses. As in Paper I, the cooling model implies 
rather strong proton SF and weak ■^P2 neutron SF in the 
NS core. In addition, we have included the effects of ^Sq 
neutron SF in the crust and found that it does not affect 
noticeably interpretation of IE 1207-52, RX J0002+62, 
Vela, PSR 0656-fl4, and Geminga as NSs with masses 
M > Md. However, this SF is crucial for interpretation of 
RX J0822-43, PSR 1055-52, and RX J1856-3754 (hotter 
for their ages) as NSs with M < Md. 

The interpretation requires the density profiles of 
Tcnsip) with rather flat and not too low maximum (<; 
5 X 10^ K) and steep decrease of Tcns(p) in the wings. 
These requirements constrain the models of nn interac- 
tion in dense matter. We have compared the ^Sq SF gaps 
A(/cFn) = ?cns/0.57 for our models 1-3 with the results 
of numerous microscopic calculations, shown in Fig. 7 in 
Lombardo & Schulze (2001). The gaps for models 1 and 
3 are typical for microscopic theories in which medium 
effects reduce ^Sq pairing not too strongly. The gap for 
model 2 is less ordinary (has too sharp decreasing slope 
at large A:Fn)- 

Note, that our results depend crucially on the choice of 
Tcp(p) at p ~ pd- By taking Tcp(p) shown in Fig. 2, we ob- 
tain the NS masses in the range from 1.1 to 1.47M0, 
in a good agreement with the well-known masses o f radi o 
pulsars in binary systems (Thorsett & Chakrabarty 1999|) . 
Adopting another Tcp{p) profile with smoother decrease of 
Tcp at high density, we would obtain higher masses for the 
same sources. The inferred mass range is also affected by 



the choice of EOS in the NS core (Paper I). However, all 
these features do not affect our conclusions on the prop- 
erties of ^Sq neutron SF needed for the NS models with 
M < Mb. 

One can see that the cooling of NSs with M < Md is 
sensitive to the density profile of free neutrons near the 
crust bottom and neutron drip point. We have used only 
one model of free-neutron distribution in the crust, assum- 
ing atomic nuclei to be spherical.lt would be interesting 
to consider the models with non-spherical nuclei at the 
crust bottom (e.g., Pethick & Ravenhall 1995) and take 



into account SF of nucleons confined in atomic nuclei. 

Let us stress that inferring from observational data 
is the most co mplica ted problem (as discussed partly in 
Yakovlev et al. 199£ ). It would be important to improve 
the limits on and the NS ages in the future high- 
quality observations; they may change considerably. Our 
interpretation is most sensitive to the data on RX J0822- 
43, and especially on PSR 1055-52, and RX J1856-3754. 
These sources seem to be excellent for testing the cooling 
theories. 
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